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ABSTRACT 

Although very successful in explaining the observed conspiracy between the baryonic 
distribution and the gravitational field in spiral galaxies without resorting to dark 
matter (DM), the modified Newtonian dynamics (MOND) paradigm still requires DM 
in X-ray bright systems. Here, to get a handle on the distribution and importance of 
this DM, and thus on its possible form, we deconstruct the mass profiles of 26 X-ray 
emitting systems in MOND, with temperatures ranging from 0.5 to 9 keV. Initially 
we compute the MOND dynamical mass as a function of radius, then subtract the 
known gas mass along with a component of galaxies which includes the cD galaxy 
with M/Lk = 1. Next we test the compatibility of the required DM with ordinary 
massive neutrinos at the experimental limit of detection [m v = 2 eV), with density 
given by the Tremaine-Gunn limit. Even by considering that the neutrino density 
stays constant and maximal within the central 100 or 150 kpc (which is the absolute 
upper limit of a possible neutrino contribution there), we show that these neutrinos 
can never account for the required DM within this region. The natural corollary of 
this finding is that, whereas clusters (T> 3 keV) might have most of their mass 
accounted for if ordinary neutrinos have a 2 eV mass, groups (T< 2 keV) cannot 
be explained by a 2 eV neutrino contribution. This means that, for instance, cluster 
baryonic dark matter (CBDM, Milgrom 2007) or even sterile neutrinos would present 
a more satisfactory solution to the problem of missing mass in MOND X-ray emitting 
systems. 
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1 INTRODUCTION 

Recently obtained data from the cosmic microwave back- 
ground (e.g., Spergel et al. 2007), large scale structure and 
type la supernovae would have us believe we inhabit a 
Universe dominated by some hitherto experimentally un- 
detected dark matter requiring an extension to the current 
standard model of particle physics (cold dark matter, CDM) 
and something called dark energy, which no one can claim 
to understand, and which drives the acceleration of the ex- 
pansion of the Universe (e.g., Perlmutter et al. 1999; Carroll 
2004; Chernin et al. 2007). However, observations appear 
to thwart the extrapolation of this ACDM theory from its 
cosmological basis to the galaxy scale because they contra- 
dict an ever growing list of CDM predictions (e.g., Moore 
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et al. 1999; Gentile et al. 2004; Kuzio de Naray et al. 2006; 
Bournaud et al. 2007). The latest gainsay of ACDM showed 
that not only can the predicted DM halos not be extrap- 
olated down to galaxy scales (where myriad clauses in the 
unknown feedback of AGNs and supernovae provide free- 
dom) but they conflict with massive clusters, where strong 
lensing shows the DM to be more concentrated than simu- 
lations predict (Broadhurst & Barkana 2008). However, DM 
concentration- virial mass relations obtained from X-ray ob- 
servations of relaxed groups and clusters agrees with ACDM 
simulations (Buote et al. 2007). 

Actually, the observed conspiracy between the mass 
profiles of baryonic matter and dark matter at all radii in 
spiral galaxies (e.g. McGaugh et al. 2007; Famaey et al. 
2007a) is more supportive of modified Newtonian dynamics 
(MOND, Milgrom 1983abc), a paradigm postulating that 
for accelerations below a ~ 10 -10 ms -2 the effective grav- 
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itational attraction approaches (gNO.o) 1 '" 2 where is the 
usual Newtonian gravitational field. 

Without resorting to CDM, this simple prescription is 
known to reproduce galaxy scaling relations such as the 
Tully-Fisher (e.g., McGaugh 2005) and Faber- Jackson (e.g., 
Sanders 2008), and the fundamental plane. More generally, it 
naturally explains the observed aforementioned conspiracy 
between the distribution of baryons and the gravitational 
field, in both low-surface brightness and high-surface bright- 
ness spiral galaxies. For reviews of serious work on MOND, 
see Sanders & McGaugh (2002), Bekenstein (2006), and Mil- 
grom (2008). 

Recently, MOND has been pushed well beyond its orig- 
inal borders with the rotation curves of several tidal dwarf 
galaxies observed by Bournaud et al. (2007), requiring new 
molecular disk DM in Newtonian gravity but conforming 
strictly to the predictions of MOND (Gentile et al. 2007 and 
Milgrom 2007a). Furthermore, taking into account the pos- 
sibility of varying the velocity dispersion anisotropy, MOND 
was shown to agree remarkably well with the line of sight ve- 
locity dispersions of the SDSS satellites (Angus et al. 2008), 
and the eight nearby MW dwarfs (Milgrom 1995; Angus 
2008, but see the individual papers for caveats). Moreover, 
the theory successfullly predicted the local galactic escape 
speed from the solar neighbourhood (Famaey, Bruneton & 
Zhao 2007b; Wu et al. 2007), and the first realistic simula- 
tions of galaxy merging were carried out (Nipoti et al. 2007; 
Tiret & Combes 2007). 

Fundamental developments in the theory of gravity have 
also added plausibility to the case for modification of gravity 
through the work of, e.g., Bekenstein (2004), Sanders (2005), 
Zlosnik et al. (2006, 2007a), Zhao (2007), and Skordis (2008), 
who have all presented Lorentz-covariant theories of grav- 
ity yielding a MOND behaviour in the appropriate limit by 
means of a dynamical vector field, which might arise from di- 
mensional reduction of a higher-dimensional gravity theory 
(Bekenstein 2006; Mavromatos & Sakellariadou 2007), or 
more generally from the fact that quantum gravity could de- 
fine a preferred rest frame at the microscopic level (Jacobson 
2008). Although rather fine-tuned and being somewhat of a 
long shot at producing a fundamental theory underpinning 
the MOND paradigm (see e.g. Bruneton & Esposito-Farese 
2007), these theories allow for new predictions, especially 
regarding cosmology. Surprisingly, such a vector field has 
precisely been shown to generate the instability that may 
produce large cosmic structures today in a MOND Universe 
(Dodelson & Liguori 2006; Skordis et al. 2006; Halle & Zhao 
2007; Zlosnik et al. 2007b) . However, before trying to build a 
consistent cosmology (see Milgrom 2008 for possible alterna- 
tive roots of MOND into cosmology) , the MOND paradigm 
has yet to completely solve the mass discrepancy in X-ray 
emitting systems. 

Groups of galaxies were studied by Milgrom (1998, 
2002) by checking the stellar mass to light ratio required 
for consistency with the line of sight velocity dispersion of 
the group galaxies w.r.t. the centre of mass. This method 
found mass to light ratios of around a solar unit which bears 
the hallmark of no mass discrepancy, however, it could not 
probe the need for dark matter at smaller radii than these 
presumably largely separated galaxies. Of course, it is most 
probable that this absence of discrepancy is correct in the 
context of low mass groups with no detected X-ray emission, 



but in the case of X-ray bright groups and clusters a better 
gauge of dynamical mass is found using measurements of 
the properties of the X-ray gas (e.g., Gerbal et al. 1992; The 
& White 1998; Sanders 1994, 1999, 2003, 2007; Aguirre et 
al. 2001; Pointecouteau & Silk 2005), or the combination of 
weak and strong lensing (e.g., Angus et al. 2007; Takahashi 
& Chiba 2007; Famaey et al. 2007c; Milgrom & Sanders 
2008; Ferreras et al. 2008). The result is that MOND can- 
not as yet explain the mass discrepancy in X-ray emitting 
clusters of galaxies, and especially in their cores, while the 
gravitational lensing map of the bullet cluster has provided 
an extremely important constraint on the nature of all the 
missing mass, i.e. that it must be of a collisionless nature 
(Clowe et al. 2006; Bradac et al. 2006; Angus et al. 2007). 

It has been conjectured (Sanders 2003, hereafter S03) 
that the mass discrepancy in galaxy clusters might be re- 
solved by the addition of a component of massive neutrinos 
m„ ~2eV, very near their maximum experimentally derived 
limiting mass of 2.2eV. These were indeed shown to be po- 
tentially consistent with the majority of clusters with tem- 
perature greater than 4 keV by S03, and with the bullet clus- 
ter (Angus et al. 2007). This hypothesis has the great advan- 
tage of naturally reproducing the proportionality of the elec- 
tron density in the cores of clusters to T 3 / 2 , as well as global 
scaling relations (Sanders 2007, hereafter S07). However, in a 
recent survey, Pointecouteau & Silk (2005, hereafter PS05) 
studied a large sample of hot clusters (>4keV) and found 
that the central density of the dark matter was generally 
greater than allowed by the Tremaine-Gunn limit on neu- 
trino density. This result was not damning, though, since 
the dynamical mass that could be accounted for by neutri- 
nos under the Tremaine-Gunn limit was generally more than 
90%. 

On the other hand, something that has never been 
addressed in the literature is the application of MOND 
to X-ray emitting groups and cool clusters in the range 
0.5< T <3.0keV. The closest anyone has come is the study 
of the large elliptical NGC 720 by Buote & Canizares (1994, 
2002) which has T~0.6keV. They showed that the major 
axis of the isophotes of the X-ray emission were signifi- 
cantly offset from the major axis of the galaxy requiring 
dark matter of at least 4 times the mass of the visible mat- 
ter. NGC 720 is currently the only elliptical galaxy for which 
this misalignment can be attributed to dark matter with rea- 
sonable confidence. A program to search for other very iso- 
lated, flattened elliptical galaxies that are sufficiently bright 
in X-rays for this type of study is underway. 

Here we primarily use the high quality galaxy group 
data from Gastaldello et al. (2007), while the data for 
A2589 and NGC 4125 come from Zappacosta et al. (2006) 
and Humphrey et al. (2006) respectively. Our data for 
higher temperature clusters come from the published work 
of Vikhlinin et al. (2006). With these 26 systems, we per- 
form a systematic check of the mass profiles of the MOND 
dark matter required for hydrostatic equilibrium of the X- 
ray emitting plasma, by subtracting the X-ray gas and the 
brightest cluster galaxy (BCG). We explore the correlation 
between the mass discrepancy in MOND and the masses and 
temperatures of the studied groups and clusters. We present 
the data in §2 and derive the MOND dynamical and dark 
masses in §3. In §4, we substract the maximum contribu- 
tion of conjectured 2eV neutrinos, by considering that the 
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neutrino density stays constant, which is bound to overesti- 
mate the contribution that neutrinos might make. Even so, 
we confirm in §5 that the central cores of hot clusters have 
densities that are too large to be consistent with ordinary 
2eV neutrinos (PS05). More significantly and in support of 
the findings of Buote & Canizares (1994, 2002), we show 
that all cool groups (0.5< T <3.0keV) have their MOND 
dynamical mass completely unaffected by these 2eV neutri- 
nos, because groups are too cool to allow them to cluster 
densely enough even at the outskirts (~ 200 kpc). In §5.3, 
we discuss other possible solutions to solve the MOND mass 
discrepancy problem, notably cluster baryonic DM (CBDM, 
Milgrom 2007b, hereafter M07) and light sterile neutrinos in 
§5.3 (e.g., Gentile, Zhao & Famaey 2008). 



2 DATA 

We use results obtained for relaxed clusters primarily from 
two recent studies of X-ray clusters. These systems were 
chosen to be among the most relaxed known in order to 
ensure that hydrostatic equilibrium holds most accurately. 

We use the published results of Vikhlinin et al. (2006) 
for 8 of the most massive clusters (2keV < T < 9keV). 
For low-mass clusters (groups) we use results for 16 groups 
of galaxies (1-2 keV) from the recent work of Gastaldello 
et al. (2007) and another low-mass cluster (A2589) from 
Zappacosta et al. (2006). Finally, we include the elliptical 
galaxy NGC 4125 (T ~ 0.5 keV) from Humphrey et al. 
(2006) to extend our sample of X-ray emitting systems down 
to the galaxy scale. 

The objects in our sample are listed in Table 1. We refer 
the reader to the references for details on the contruction of 
the density, temperature, and Newtonian mass profiles. We 
note that for the massive systems in the Vikhlinin et al. 
sample we exclude the central ~ 20 kpc from the analysis. 
We do this for consistency with their study, even though 
some of their systems (most notably A2029), do not exhibit 
substantial morphological irregularities in their cores that 
would lessen the validity of the approximation of hydrostatic 
equilibrium. 



3 MOND DYNAMICAL MASS IN GROUPS 
AND CLUSTERS 

To compute gravitating masses we assume the intracluster 
medium is represented by a spherical single-phase ideal gas 
in hydrostatic equilibrium. Although the X-ray isophotes 
of relaxed clusters are approximately circular with mod- 
est ellipticity, the underlying mass distribution is inferred 
to have substantial ellipiticy (0.4-0.6; Buote & Canizares 
1996). Nevertheless, many previous studies have shown that 
assuming spherical symmetry for relaxed clusters introduces 
fairly small errors < 20% which are acceptable for our pur- 
poses (e.g., Tsai, Katz, & Bertschinger 1994; Navarro et al. 
1995; Buote & Canizares 1996; Evrard et al. 1996; Gavazzi 
2005). 

We derive the MOND dynamical mass for each system 
following the approach of Sanders (1999). From the tem- 
perature and density profiles of the hot gas, the centripetal 



gravitational acceleration, g, can be deduced from the equa- 
tion of hydrostatic equilibrium independently of the gravi- 
tational theory: 



g(r) = 



-fcT(r) 



dlnp x (r) d In T(r) 
dlnr 



dlnr 



(1) 



where w is the mean molecular weight and wm p = 5.2 x 
10 -58 Mq and the combination kT(r) is in units of keV. 
From this, the Newtonian dynamical mass can easily be de- 
duced. The studies cited in the previous section, from which 
we obtain our data, use an NFW profile (Navarro, Frenk & 
White 1995) to match the dark matter profile. 

Note that since we model clusters as spherical systems, 
we can ignore the curl field of MOND (Bekenstein & Mil- 
grom 1984; Ciotti, Londrillo & Nipoti 2006; Angus, Famaey 
& Zhao 2006), and since the external field effect from Large 
Scale Structure (e.g., Angus & McGaugh 2007; Famaey et 
al. 2007b; Wu et al. 2007) is much smaller than the typ- 
ical gravitational acceleration in the region of interest, we 
can simply use the relation of Milgrom (1983a) to derive the 
MOND dynamical mass: 



M m {r) = g5(r)M[sW/a„] 



(2) 



where G = 4.42 x 10 _3 pc(kms" 1 ) 2 Mg 1 is Newton's grav- 
itational constant. Here the function [i{x) is chosen to be 
the simple function shown to have a good fit to the terminal 
velocity curve of the Milky Way (Famaey & Binney 2005) 
and to rotation curves of external galaxies (Famaey et al. 
2007a; Sanders & Noordermeer 2007): 



H(x) = x/(l + x), 



(3) 



while taking fj, — 1 just reduces to the Newtonian dynamical 
mass. This means that the MOND dynamical mass M m is 
related to the Newtonian one M n by 



M„ 



M n (r) 



l + a /g(r)' 



(4) 



so that in MOND, there is a truncation of the dynam- 
ical mass M m (r) at low accelerations. Here we take the 
MOND acceleration constant to be a = 3.6 (kms _1 ) 2 /pc = 
1.2 x 10~ 10 ms -2 . Moreover, from Eq. ((4}, one can see that at 
large radii, when d\n M n / d\nr < 1 (as is the case for NFW 
profiles) and g(r) « a , then the MONDian enclosed dy- 
namical mass begins to drop with radius. This is of course 
unphysical, and is caused by MOND generically predicting 
a logarithmic potential at large radii. This makes studies 
that do not consider the entire dynamical mass profile at 
all radii worthless when based on a Newtonian NFW pro- 
file. In this work, we circumvent this problem by flattening 
the enclosed MOND mass profile at its maximum value, and 
numbers in Table 1 are quoted at the radius corresponding 
to this maximum value. Further away, MOND predicts that 
the potential should not conform to an NFW profile but 
rather to an isothermal sphere. 

For the objects in our sample from Gastaldello et al. 
(2007), Humphrey et al. (2006), and Zappacosta et al. 
(2006), we take the Newtonian gravitating mass profiles used 
in those studies and simply extract the MOND dynami- 
cal mass using Eq. ((4]). Those studie produced 20-30 Monte 
Carlo simulations of Newtonian mass profiles for each object, 
which were used to produce the errors. For MOND profiles, 
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because the /i-function is not precisely known, any errors 
on the data are dwarved by the intrinsic errors. Therefore, 
we take the error on the MOND mass as half the difference 
between M m calculated with the simple /i-function (Eq. [3] 
which gives a relatively low mass) and M m with the stan- 
dard /i- function (ji = x/\/l + x 2 , which gives a relatively 
high mass). These limits contain the majority of the reason- 
able MOND profiles and so can be naively taken as 1-a. 

Since we do not have error estimates for the density 
and temperature profiles of the objects in the Vikhlinin et 
al. sample, we just use the best-fitting values. However, as 
the objects are more massive and generally brighter than 
the lower mass objects, we expect the relative uncertainty 
on the massive systems to be typically comparable to or less 
than the lower mass systems. In any case, the systematic 
errors on the MOND masses calculated by using different \i- 
functions are large enough to account for any small random 
errors in the Newtonian mass. 



3.1 Subtracting the X-ray gas and galaxies 

We integrated the plasma component to find the observed 
enclosed gas mass, Mx{r) from 



M x (r) 



4:irpx(r)r dr 



(5) 



where px(r) is the density of plasma. 

The empirical scheme for the galaxy mass pro- 
file employed by S03 and PS05, i.e. M*(r)/M x (r) « 
0.4(fcT(r)/keV) _1 completely neglects the contribution of 
galaxies to the cluster mass below f 50 kpc where the most 
massive galaxies reside. Furthermore, for groups this is a 
poorly motivated scheme and so for them we only include 
the mass of the BCG. This is reasonable as the MOND dy- 
namical mass of groups saturates at around 100 kpc and we 
checked the surrounding 100 kpc of the BCG and found no 
galaxies contributing more than 5% which is easily contained 
within the errors on M/Lk- 

For the clusters we used both the BCG luminosity and 
the empirical galaxy density of S03 and PS05. The luminos- 
ity of the BCG in the K band from the Gastaldello et al. 
(2007) sample are taken from their paper, while those from 
Vikhlinin et al. (2006) are taken from Lin & Mohr (2004) 
and the 2MASS survey. All are employed using a Hernquist 
profile. We subtracted this using M/Lk~1 which is stan- 
dard for these very old stellar populations and a Kroupa 
IMF (Kroupa 2001); however, if a Salpeter IMF is used, we 
can expect M / Lk to increase by around 50% (Humphrey 
2006). The X-ray studies of groups and clusters (Gastaldello 
et al. 2007; Zappacosta et al. 2006) typically find M/L K <1, 
sometimes << 1 as they try to incorporate an NFW halo. 
By using M/Lk = 1, we are making a conservative estimate 
of the amount of missing mass. 

Subtracting both the X-ray mass and the galaxies from 
the total MOND dynamical mass leaves the total dark mass 
(Mdm) which is simply the mass of dark material necessary 
to reach agreement with the MOND dynamical mass. 



4 ORDINARY NEUTRINOS AS THE MOND 
DARK MASS? 

At least two of the three active neutrinos (v e , v T and v^) 
have non-zero masses, meaning that they must be part of 
the mass budget of the Universe. Interestingly, in ACDM 
cosmology, it is possible to put stringent limits on the sum of 
the three neutrino masses, from the angular power spectrum 
of the Cosmic Microwave Background and from the slope of 
the matter power spectrum (e.g., Zunckel & Ferreira 2007; 
Host et al. 2007), essentially because a neutrino component 
that is too massive would not leave enough room for CDM 
in the matter budget of the Universe, resulting in the loss 
of small-scale power. 

However, these constraints are not necessarily valid in 
a modified gravity cosmology (e.g., McGaugh 2004; Skordis 
et al. 2006) where the stronger gravity (covariantly in the 
form of an additional vector field) can effectively play the 
role of cosmic dark matter (thus not necessarily tracking the 
baryons, contrarily to the non-covariant case) by adding a 
new instability in the process of structure formation (e.g., 
Bekenstein 2004; Dodelson & Liguori 2006; Halle & Zhao 
2007; Zlosnik et al. 2007b). 

Alternatively, model-independent experimental limits 
on the electron neutrino mass from the Mainz/Troitsk ex- 
periments, counting the highest energy /3-decay electrons 
of 3 H -» 3 He + +e~ + 18.57 keV (the more massive 

the neutrinos, the lower the cutoff energy of electrons), are 
m„ < 2.2eV\ The KATRIN experiment (under construc- 
tion) will be able to falsify 2eV neutrinos at 95% confidence 
within months of taking data in 2009. 

If the neutrino mass is substantially larger than the 
mass differences, then all types have about the same mass, 
and the cosmological density of three left-handed neutrinos 
and their antiparticles (e.g., S07) would be 



H„ = 0.062m„ 



(6) 



where m v is the mass of a single neutrino type in eV. If one 
assumes that clusters of galaxies respect the baryon-neutrino 
cosmological ratio, and that the MOND DM is mostly made 
of neutrinos as suggested by S03 and S07, then the mass of 
neutrinos must indeed be around 2 eV. 

In their modelling of the CMB anisotropies, Skordis et 
al. (2006) showed that such a component of 2eV neutrinos 
could actually prevent the MOND Universe from acceler- 
ating too much, and could thus yield the correct angular- 
distance relation to get the correct position of the peaks in 
the angular power spectrum of the CMB. 

The main limit on the neutrino ability to condense in 
clusters comes from the Tremaine-Gunn limit (Tremaine & 
Gunn 1979), stating that the phase space density must be 
preserved during collapse. This is a density level half the 
quantum mechanical degeneracy level. The maximum den- 
sity for a cluster of a given temperature, T, is defined for a 
given mass of one neutrino type as 



7 x l0- 5 M Q pc- 



(— V 

\lkeVJ 



f m v y 
\2eVJ 



(7) 



S03 showed that such 2 eV neutrinos at the limit of detection 
could indeed account for the bulk of the dynamical mass in 
his sample of galaxy clusters of T > 4keV (see his Fig. 8). 
Angus et al. (2007) also showed that such neutrinos could ac- 
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Figure 1. Shows the mass profiles of the components of the dynamical mass for 4 representative clusters: N5044 (T=1.0 keV), N533 
(T=1.2 keV), A2717 (T=2.2 keV), and A2029 (T=8.5 keV). The total MOND dynamical mass (M,„) is in black with thick-solid linetype 
and its 1-cr error is the thin-solid black line (there is virtually no error below). The Newtonian dynamical mass is black with a dotted 
linetype. The solid red line corresponds to the observed mass of X-rays and the dotted red line is the mass of the BCG if M/Lx=l for 
which we use a Hernquist profile. The green line is the maximum necessary contribution of neutrinos; by this we mean it is maximal 
(Eq[7J in the centre where the dynamical mass is unexplained, whereas at the outskirts the neutrinos no longer have maximum density, 
instead they have the density necessary to complete the budget after gas and the residual mass at the centre has been accounted for. The 
thick solid blue line is the residual mass unexplained by the neutrinos and gas and its 1-<t error is the thin solid blue line. The dotted 
blue line is the necessary residual mass if we have no significant neutrino density and again the 1-cr error is the thin dotted blue line. Of 
course there is another error associated with the mass of the neutino, but we have fixed the neutrino mass at 2eV. Obviously a range 
of blue solid lines are possible between the solid blue and dotted blue depending on neutrino mass. If the neutrino has a mass <<2eV 
then we simply recover the dotted blue line. Clearly, one sees that for low temperature clusters such as N5044, no DM is present at r 
> 150 kpc. The errors on the different lines are not independent from each other. The firm result here is that no random or systematic 
errors can ever make the green line (which overestimates the true contribution of 2 eV neutrinos by assuming a fully constant density in 
the core) reach the thick dotted blue line in the core of clusters, or even anywhere in groups, if the mass of ordinary neutrinos is smaller 
than its experimental upper bound. 
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count for the weak lensing map of the bullet cluster (Clowe et 
al. 2006; Bradac et al. 2006), while S07 maintained that this 
hypothesis has the great advantage of naturally reproduc- 
ing the proportionality of the electron density in the cores of 
clusters to T 3 ^ 2 , as well as global scaling relations. However, 
looking at the central region of clusters, PS05 showed that 
neutrinos could not account for the dark matter all the way 
to the centre because the density reaches values much larger 
than the Tremaine-Gunn limit. However, this residual mass 
is only a few percent of the dynamical mass explainable by 
neutrinos. 

Here, we choose as a conservative approach to take the 
temperature of the neutrino fluid as being equal (due to vio- 
lent relaxation) to the mean emission weighted temperature 
of the gas. We also assume that they contribute maximally, 
i.e. that their density is given by the Tremaine-Gunn limit, 
and constant (which is obviously untrue since the fluid obeys 
the equation of state of a partially degenerate neutrino gas, 
see e.g. S07). We thus overestimate the true contribution 
that neutrinos might make. This is of little concern, how- 
ever, since, as we shall see in the next section, we do not wish 
to conclude anything beyond the fact that even overestimat- 
ing the physical contribution of 2 eV neutrinos to the mass 
budget of the cluster does leave a massive, dense central 
component of unexplained mass, and contributes nothing to 
groups which again require a dense residual component. We 
shall estimate the importance of this residual component in 
clusters and groups, but beyond that we shall also consider 
the amount of DM needed if neutrinos have a negligible mass 
and do not contribute at all to the mass budget. 



5 RESIDUAL MASS 

In Fig. [1] we plot for four representative objects several com- 
ponents of the group or cluster mass distinctly as functions 
of radius. The total MOND dynamical mass (M m ) is in black 
with solid linetype, the Newtonian dynamical mass (M n ) is 
black with a dotted linetype. The red line corresponds to 
the observed mass of X-rays (Mx) and the dotted red line 
is the mass of the BCG assuming M/Lk=1 and a Hernquist 
profile. The green line is the maximum necessary contribu- 
tion of neutrinos (M„). The solid blue line is what we call 
the residual mass (M m -x-v-*) unexplained by the neutri- 
nos, galaxies and gas, whereas the dotted blue line is the 
necessary total dark mass (M m -x-* = Mom) if we have 
no significant neutrino density. The various thin lines cor- 
respond to the systematic errors coming from the choice of 
/^-function. The four representative clusters have respective 
temperatures 1.0, 1.2, 2.2 and 8.5 keV. Clearly, one sees 
that for low-temperature groups such as NGC5044, no DM 
is present at r > 150 kpc. 

We conclude that neutrinos might explain the mass dis- 
crepancy in the outer parts of the groups/clusters, especially 
those with high temperatures (T > 3 keV), keeping in mind 
that no dark mass at all is necessary in the outer parts for 
the cooler clusters. However, to examine the real contribu- 
tion of massive neutrinos in the outer parts of rich clusters, 
one should describe the equilibrium distribution expected 
for the neutrinos, and then examine whether this physically 
well-motivated profile can fit the data. This is however a 
non-trivial task in the context of MOND, needing to solve a 



Lame-Emden-like equation for polytropic models. Moreover, 
if the models are made more realistic by not considering 
purely self-gravitating neutrino spheres but also including a 
component of baryons, the task becomes even less trivial. 
This is far beyond the scope of the present paper, and will 
be the subject of further studies. 

Nevertheless, the important result is that even with our 
overestimate of the neutrino contribution, we found that 
neutrinos cannot account for the dark mass in the central 
150 kpc. Much fuss has been made in the literature about the 
neutrinos being unable to account for this core of residual 
mass (PS05; Takahashi & Chiba 2007, Ferreras et al. 2008), 
which is true. Indeed, every one of the sample clusters has 
a dominant central core of residual mass. Nevertheless, al- 
though all these systems have a residual dark component, 
it has been neglected as to how serious this extra compo- 
nent is, which the neutrinos cannot account for, and how 
this varies with temperature and mass of the cluster. Actu- 
ally, in Fig. [4] which plots the ratio of residual mass (with 
and without neutrinos) to MOND dynamical mass for all 
clusters, we see that for the hot clusters studied by S03 and 
PS05 the residual mass (after neutrinos have been added) 
is a small fraction, with gas and neutrinos being more im- 
portant globally (even though, as stated above, the contri- 
bution of neutrinos is a bit overestimated here) . Conversely, 
for groups like NGC5044, the residual mass is completely 
dominant and is up to four times more significant than the 
stellar and gas components. Recall that the neutrinos can- 
not lay claim to the high densities in the cores due to the 
Tremaine-Gunn limit (and low temperatures of groups), but 
at a certain radius begin to dominate the dynamical mass. 
Given that our approximation is bound to overestimate the 
neutrino contribution, this radius is actually the minimal ra- 
dius at which neutrinos might begin to dominate the mass 
budget of the cluster. This radius r u is plotted as a function 
of temperature in Fig[2] for 2eV neutrinos. It is intriguing 
that if we do live in a universe filled with ~2eV neutrinos 
then the core sizes of residual DM in large and small clus- 
ters are very similar. Of course, if neutrinos are not present, 
then the residual DM is by far the most dominant compo- 
nent, moreso than the gas and galaxies for all groups and 
clusters. It is also important to realise that any departures 
from hydrostatic equilibrium generally cause underestimates 
of the total mass; i.e., from non-thermal pressure support. 
This would only exacerbate the problems highlighted in this 
paper. 

In Table 1 we list the different mass components in the 
sample of clusters, and in Fig[3]we plot the MOND dynami- 
cal mass of the clusters against temperature, which is fitted 
with a curve of the form M m oc T 2 ' 4 , close to the MOND 
prediction (S07) of a relation M m oc T 2 . In order to get 
a strict M m oc T 2 the system has to reach the asymptotic 
MOND regime and the T used should be the asymptotic 
T, not the emission weighted average (or otherwise) and of 
course, must be constant. Furthermore, the systems should 
have the same logarithmic density slope at the measured 
radius. 

We show in Fig(4]the fraction of residual mass as a func- 
tion of cluster mass. If neutrinos are not present, then the 
fraction of residual dark mass is constantly high (between 
60% and 80%) because it must account for the largest chunk 
of the dynamical mass, however, if neutrinos are present 
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Figure 2. Shows the radius below which neutrinos do not con- 
tribute in sufficient density to the MOND hidden mass. 



then the residual mass becomes less important for heavier 
clusters. In Fig[5] one sees that the total dark mass scales 
with temperature, but that if 2 eV neutrinos are present 
the residual mass appears to saturate at a temperature of 
around 3 keV and a mass of ~ 10 13 M©. Interestingly, at 
the opposite end of the scale, the residual mass is steeply 
falling towards zero for groups of T < 2keV, meaning that 
the amount of residual mass could in some sense be linked 
with the energy of X-rays themselves. Below the 0.5 keV 
threshold we get back to systems like the groups studied by 
Milgrom (1998, 2002), where no dark mass is needed and no 
X-rays are detected. 



5.1 Brightest Cluster Galaxies 

We express the quantity of residual mass in the central parts 
of the cluster as a MOND dynamical mass-to-light ratio for 
the BCG of the cluster, ranging between 1.7 and 20 (see 
Fig[7|. This mass-to-light ratio corresponds to an amount 
of mass that must be present within the radius r„. This 
high dynamical mass-to-light ratio might perhaps be due 
to the interaction of this BCG with the neutrino fluid, al- 
lowing the fluid to heat and the neutrinos to clump more 
densely. From Fig(5]we see that the residual mass does in- 
deed scale with the K-band luminosity of the BCG which is 
as expected if the neutrino fluid is heated allowing higher 
densities. If it was the case that the BCGs had large M/Lk 
we should be able to fit the stellar mass profile (with variable 
normalisation) without a DM profile and it should describe 
the data well. Unfortunately, the mass profiles of galaxy 
groups/clusters are not well fit by a De Vaucoulers profile 
with R e set to that of the central galaxy but allowing for 
variable M/Lk- If we did this, we would over predict the 
interior mass profile. In addition, the barycentres of many 
clusters do not coincide with the BCG. This scaling of the 



Figure 3. Shows the scaling of the MOND dynamical mass with 
cluster temperature. The fitted line is M m cx T 24 . This is higher 
than the MOND prediction of M m oc T 2 for the reasons explained 
in (0 



residual mass with the luminosity of the BCG could thus 
rather mean that the heaviest BCGs reside in clusters with 
more X-ray emitting hot gas, and that large amounts of cold 
gas clouds are in turn also present (see § 6.2). 



5.2 Transition temperatures 

There is an increase in significance of the residual mass com- 
ponent left after the subtraction of gas and neutrinos for 
lower temperature clusters. Furthermore, for groups below 
2keV, the contribution of neutrinos is negligible compared 
to the X-ray gas. Conversely, neutrinos have the potential 
to be of great importance in the dynamics of clusters hot- 
ter than 2keV, and subsequently are likely to help seed the 
collapse of these structures from cosmological perturbations 
(see Sanders 2008 for an introduction to galaxy formation 
in MOND). 

It is well known that Romanowsky et al. (2003) 's stellar 
dynamical study of stellar tracers using planetary nebulae 
in purposely X-ray faint elliptical galaxies, generally do not 
require dark matter in the central parts (which is known 
as the problem of the dearth of dark matter in elliptical 
galaxies). Consequently, it is not surprising that they have 
also been shown to be consisent with MOND (Milgrom & 
Sanders 2003, Tiret et al. 2007). In addition, the groups of 
galaxies studied by Milgrom (1998,2002) appear consistent 
with no dark matter. However, those studies suffer from the 
well-known degeneracy from slight freedom in the M/L and 
velocity dispersion anisotropy, which allows a wide range of 
mass profiles to be consisent with the data. 

Angus et al. (2008) used the two cleverly assembled 
mock galaxy groups of Klypin & Prada (2007) and computed 
the necessary M/L of the central galaxy in producing the los 
velocity dispersions of the stacked satellites. The lower mass 
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Figure 4. Shows the fraction of dark mass (red) and of residual 
mass after neutrinos have been taken into account (black) as a 
function of the MOND dynamical mass. Clearly the ratio of dark 
mass stays constantly high (>0.6) if neutrinos are not contribut- 
ing, but the fraction of residual mass becomes less important for 
heavier (and hotter) clusters when neutrinos are added. 



Figure 6. Shows the marginal correlation of the K-band lumi- 
nosity of the BCG with cluster temperature. 
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Figure 5. Shows the residual mass discrepancy vs. temperature 
after subtraction of X-ray gas, the BCG and of the maximum 
contribution of neutrinos (black) and for the case with no neu- 
trinos (red). In the case of significant neutrino contribution, this 
discrepancy is a fairly constant value of ~ 1O 13 M0 for T > 3 keV 
but drops steadily to zero for T < 3 keV. If neutrinos are not 
present, the dark mass continues to rise with temperature. 




Figure 7. Shows the necessary K-band M/L ratio of the BCG 
to explain the residual central mass discrepancy when neutrinos 
are present. It ranges between 2 and 20, meaning the stellar pop- 
ulation of the BCG is probably not enough. If no neutrinos are 
present, then the ratio > 100 for some clusters (see the table). 



host is around 8x10 Mq using a sensible M/L in the g- 
band of 1.9-3.3, whereas the higher luminosity host requires 
2.8xlO n M which can only be achieved with a M/L of 4.2- 
6.6 whereas the acceptable range is 3-5. So, there exists some 
evidence that the higher mass elliptical requires a moderate 
amount of DM. 

We can only be certain that HSB galaxies like those 
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Figure 8. Shows the scaling of the residual mass discrepancy 
as a function of BCG luminosity, with (black) and without (red) 
neutrinos. The scaling may reflect the deeper potential well heat- 
ing the neutrino fluid allowing a higher density far beyond that 
allowed by the mean cluster temperature, or the fact that dark 
matter in the form of cold gas is more abundant in brighter galax- 
ies. More trivially, it is just the natural expectation of brighter 
galaxies residing in more massive clusters. 



studied by Sanders & Noordermeer (2007) are fully consis- 
tent with MOND and no abundance of galactic DM. Be- 
low systems of this mass (at least) there is no need for DM 
in MOND. Studies of very X-ray luminous elliptical galax- 
ies have not generally examined the viablity of alternative 
gravity theories. A recent exception is the study of the ve- 
locity dispersions of globular clusters within the giant el- 
liptical NGC 1399 (Richtler et al. 2008), which has a two- 
temperature gas of 0.9-1.5 keV (Buote 2002), that has re- 
vealed some need for DM (within the aforementioned uncer- 
tainties on the anisotropy). Another exception is NGC 720, 
where the observed misalignment of the major axis of the 
stars and hot gas, along with the elongation of the X-ray 
isophotes, has been suggested to imply a substantial mass 
discrepancy (Buote & Canizares 1994; Buote et al. 2002). 
More recently the analysis of the large elliptical galaxy 
NGC 4125 (included here) with T=0.5keV shows more ev- 
idence for MOND DM in ellipticals. Here we find that we 
require NGC 4125 to have a minimum M/Lk=1-7 mean- 
ing, for the only system in our sample, there is less DM 
than luminous matter (note that this mass-to-light ratio in 
the K-band is obtained after subtracting the hot gas, and 
cannot be due to this component itself). 

The requirement of DM in MOND at 0.5 keV for 
NGC 4125, 0.6 keV for NGC 720 and the SDSS galaxies but 
no requirement for HSBs (Sanders & Noordermeer 2007) or 
X-ray dim ellipticals and galaxy groups (Milgrom & Sanders 
2003, Milgrom 2002), necessitates a transition temperature 
below which the fraction of residual mass diminishes in im- 
portance and above which it increases in importance, peaks 



(between 0.5 keV and 2 keV) and then falls at higher tem- 
peratures if neutrinos are present (because the residual mass 
saturates at 10 13 Mq but the dynamical mass continues to 
rise). If neutrinos do not contribute significantly, the residual 
amount of DM would directly scale with temperature. Why 
systems cooler than 0.5 keV show virtually no mass discrep- 
ancy in MOND must be related to the formation mechanism 
of such systems, and to the nature of the residual mass com- 
ponent found in X-ray groups (or of the whole MOND DM in 
groups and clusters if neutrinos have a negligible mass). Let 
us note that the very presence of DM in MOND seems to be 
synonymous with the presence of ionised gas and X-ray emis- 
sion, which might hint at the presence of similar quantities 
of unseen cold gas that could perhaps explain this residual 
mass, and maybe even the whole of the mass discrepancy 
problem of MOND (M07). 

5.3 Cold gas clouds 

In the global baryon inventory of the Universe 50% of the 
baryons produced during Big Bang nucleosynthesis (BBN) 
are still missing, and assumed to be in the warm-hot inter- 
galactic medium (WHIM, see e.g. Roncarelli et al. 2006). 
Depending on the actual amount of baryons in the WHIM, 
there is thus some freedom for the residual mass to be 
baryonic in MOND. Indeed Fukugita & Peebles (2004) es- 
timate that the observed baryons in clusters only account 
for 5% of those produced during BBN, while it is clear from 
our Fig. [1] that the gas mass is increasing faster than the 
MOND dynamical mass, and tends to a DM/baryons ratio 
of about 1 at large radii. This means that if only 45% of 
baryons are present in the WHIM, the remaining 5% would 
be sufficient to solve the whole problem of missing mass in 
MOND clusters (even if neutrinos were found by KATRIN 
in 2009 to have mass much below the present-day experi- 
mental limit). One could thus say that, while CDM puts the 
missing baryons problem in individual galaxies, MOND has 
put it in galaxy clusters (see also McGaugh 2007) . It should 
be highlighted again that DM in MOND only appears in 
systems with an abundance of ionised gas and X-ray emis- 
sion (even in the case of galaxies such as NGC 4125). It is 
then no stretch of the imagination to surmise that these gas 
rich systems have equal quantities of molecular hydrogen (or 
other molecules), in e.g. some compact form to satisfy the 
collisionless citerion imposed by the bullet cluster. 

M07 has recently proposed that the missing mass in 
MOND could entirely be "cluster baryonic dark matter" 
(CBDM), e.g. in the form of cold, dense gas clouds. There is 
an extensive literature discussing searches for cold gas in the 
cores of galaxy clusters (e.g., Donahue 2006) but what is usu- 
ally meant there is quite different from what is meant here, 
since those searches consisted in trying to find the signature 
of diffuse cold molecular gas at a temperature of ~ 30 K. 
The proposition of M07 rather relies on the work from Pfen- 
niger & Combes (1994) and Walker & Wardle (1998), where 
dense gas clouds with a temperature of only a few Kelvins 
(~ 3 K), Solar System size, and of a Jupiter mass, were con- 
sidered to be possible candidates for both galactic and ex- 
tragalactic DM. Since the CBDM considered in the context 
of MOND cannot be present in galaxies, it is however not 
subject to the galactic constraints on such gas clouds. Note 
that the total sky cover factor of such clouds in the core of 
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the clusters would be of the order of only 10 -4 , so that they 
would only occult a minor fraction of the X-rays emitted by 
the hot gas (and it would be a rather constant fraction). For 
the same reason, the chances of a given quasar having light 
absorbed by them is very small. Still, M07 notes that these 
clouds could be probed through X-ray flashes coming out of 
individual collisions between them. 

Of course, this speculative idea also raises a number of 
questions, the most serious one being how these clumps form 
and stabilize, and why they form only in clusters, groups 
and some ellipticals, but not in individual spiral galaxies. As 
noted above, the fact that missing mass in MOND is nec- 
essarily associated with an abundance of ionised gas could 
be a hint at a formation and stabilization process some- 
how linked with the presence of hot gas and X-ray emission 
themselves. Then, there is the issue to know whether the 
clouds formation would be prior or posterior to the cluster 
formation. We note that a rather late formation mechanism 
could help increase the metal abundance, solving the prob- 
lem of small-scale variations of metallicity in clusters when 
the clouds are destroyed (Morris & Fabian 2003). M07 also 
noted that these clouds could alleviate the cooling flow co- 
nundrum, because whatever destroys them (e.g. cloud-cloud 
collisions and dynamical friction between the clouds and the 
hot gas) is conducive to heating the core gas, and thus pre- 
venting it from cooling too quickly. Such a heating source 
would not be transient and would be quite isotropic, con- 
trary to AGN heating. This highly interesting and innova- 
tive scenario to alleviate the cooling flow problem should be 
the subject of further studies, but is however not the focus 
of our present paper. 

5.4 Sterile neutrinos 

In addition to the three flavours of ordinary neutrinos, there 
may be one or two sterile neutrinos of which the particle 
mass is neither constrained by the CMB or structure for- 
mation (e.g., Abazajian et al. 2001), and surely not in a 
modified gravity scenario. Such sterile neutrinos might be 
light enough (eV scale, see e.g. Maltoni & Schwetz 2007) 
not to heavily perturb MOND fits in individual galaxies, 
but heavy enough to account for the missing mass in groups 
and clusters. Gentile et al. (2008) postulated the existence 
of one or two species of 5 eV neutrinos in this context. Here, 
we computed the minimum sterile neutrino mass for a sin- 
gle species using the residual mass densities of NGC 533 and 
A2029 at 50 kpc, assuming that ordinary neutrinos have a 
negligible mass. We found that the minimal mass required 
was 5.6eV and 4.0eV respectively, meaning that the exis- 
tence of a 6 eV sterile neutrino would have the potential 
to completely solve the cluster problem in MOND. How- 
ever, further studies should describe the equilibrium distri- 
bution expected for such sterile neutrinos, and then examine 
whether physically well-motivated mass profiles can fit the 
data. 



that this effect was most important at the tidal boundary of 
merging, rapidly varying, systems, such as the high collision 
speed of the bullet cluster. It however remains to be seen 
how this effect would offer an explanation to the MOND 
mass discrepancy found for the quasi-static systems studied 
in this paper, including X-ray emitting elliptical galaxies. 
Alternatively, the introduction of an additional scalar field 
(e.g., Sanders 2005) in covariant theories of MOND can in- 
duce the presence of bosonic DM (through the scalar field 
oscilations) that cannot accumulate in spiral galaxies, but 
that could help the cause of the missing mass in groups and 
clusters. 



6 CONCLUSION 

We have used high quality data to decompose the mass pro- 
files of 26 X-ray emitting systems in MOND, with tempera- 
tures ranging from 0.5 to 9 keV. We confirmed that galaxy 
groups and clusters in MOND require a huge DM component 
to make up for all the dynamical mass (60-80% of the total 
dynamical mass at the last observed radius, dropping there- 
after). We have shown that, whatever their precise equi- 
librium distribution, hypothetical 2 eV ordinary neutrinos 
can never explain a component of residual mass enclosed 
within the central 100 or 150 kpc. This confirms the result 
of PS05, but we add an important corollary to this finding, 
namely that groups with T<J 2 keV cannot be explained by 
a 2 eV neutrino contribution. Indeed, we have probed far 
lower masses (< 1O 13 M ) than PS05, and have shown that 
this residual mass component, which is present in all known 
X-ray bright systems hotter than 0.5keV, is the dominant 
one for groups with T< 2keV. 

This leads us to believe that, if MOND is correct, the 
current baryon inventory of galaxy clusters is unfinished and 
that the unaccounted for baryons known to exist from argu- 
ments of Big Bang Nucleosynthesis (depending on the frac- 
tion residing in the WHIM) are indeed buried in the central 
regions of clusters in some hitherto undetectable form. Con- 
straints from the bullet cluster tell us that this mass must be 
compact and cool. This could hint at dense cold gas clouds 
as recently suggested by M07. It could alternatively hint at 
the existence of a 4th light sterile neutrino of ~ 6eV mass. 
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5.5 Field oscillations 

Finally, we note that it has recently been suggested that non- 
trivial effects of the vector field in covariant formulations of 
MOND might create a DM-like effect (e.g., Zhao 2007), alle- 
viating the need for neutrinos and/or CBDM. It was argued 
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Table 1. (1) Designations of the 26 groups and clusters. (2) Mean emission weighted temperature. (3) Radius of last data point. (4) 
Newtonian dynamical mass at overdensity A=1250, 2500, 200 and 500 for the superscripts a, b, c and d respectively. (5) X-ray gas 
mass at r max , may be larger than required for equilibrium. Error is always less than 1% (6) MOND dynamical mass using the simple 
/i-function . (7) Residual dark mass required if neutrinos are not present. (8) Contribution of neutrinos to M m at r max . (9) Minimal 
radius beyond which neutrinos have high enough density to supply remaining dynamical mass. (10) Residual mass if neutrinos contribute 
maximally. (11) K-band luminosity of BCG. (12) Necessary Mj to account for (10). (13) The reference from which the data is taken, 
(a) indicates Gastaldello et al. (2007a), whereas (b) is Vikhlinin ct al. (2006), (c) Zappacosta et al. (2006) and (d) Humphrey et al. 
(2006). Three clusters overlap between Vikhlinin et al. (2006) and Gastaldello ct al. (2007a) which we pair in the table and have similar 
results. 
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